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ABSTRACT 



Context. The study of the atmosphere of red supergiant stars in general and of Betelgeuse (a Orionis) in particular is of prime 
I ■ importance to understand dust formation and how mass is lost to the interstellar medium in evolved massive stars. 

Aims. A molecular shell, the MOLsphere (Tsuji, 2000a), in the atmosphere of Betelgeuse has been proposed to account for the near- 
and mid-infrared spectroscopic observations of Betelgeuse. The goal is to further test this hypothesis and to identify some of the 
CZ3 , molecules in this MOLsphere. 

Methods. We report on measurements taken with the mid-infrared two-telescope beam combiner of the VLTI, MIDI, operated between 
7.5 and 13.5 yum. The data are compared to a simple geometric model of a photosphere surrounded by a warm absorbing and emitting 
shell. Physical characteristics of the shell are derived: size, temperature and optical depth. The chemical constituents are determined 
with an analysis consistent with available infrared spectra and interferometric data. 

Results. The MIDI data are well modeled with a geometrically thin shell whose radius varies from 1.31 to 1.43 R+ across the N band 
with a typical temperature of 1550 K. We are able to account for the measured optical depth of the shell in the N band, the ISO-SWS 
spectrum and K and L band interferometric data with a shell whose inner and outer radii are given by the above range and with the 
following species and densities: H 2 (7.1 ± 4.7 x 10 19 cnT 2 ), SiO (4.0 ± 1.1 x 10 20 cnT 2 ), A1 2 3 (2.4 ± 0.5 x 10 15 cnT 2 ). 
Conclusions. These results confirm the MOLsphere model. We bring evidence for more constituents and for the presence of species 
participating in the formation of dust grains in the atmosphere of the star, i.e. well below the distance at which the dust shell is 
detected. We believe these results bring key elements to the understanding of mass loss in Betelgeuse and red supergiants in general 
f*""*. , and bring support to the dust-driven scenario. 
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d ' 1. Introduction reinterpreted previously published data (Wing & Spinrad 1970) 

and proposed a 1500 + 500 K water vapor and CO shell around 

The atmospheres of late-type supergiant stars are not a priori Betelgeuse and ^ Cepi the MOLsphere. 
supported by large amplitude pulsations of the star as is the case 

in Mira stars but are still much more extended than expected in Multi-telescope optical interferometry offers new and im- 
hydrostatic equilibrium. In the case of Betelgeuse, this is directly portant constraints on the spatial location of these regions in- 
illustrated by apparent sizes ranging from 42 up to 65 mas, the ferred from spectroscopy. |Perrin et al. (2004a)] have shown that 
size being minimum in the near-infrared and larger in the vis- the near-infrared structure of Betelgeuse can be described by two 
ible and ultraviolet. The large extent of the atmosphere is not main components: a photosphere of temperature 3641 + 53 K 
understood and makes it complex to study. and a gaseous shell at a temperature of 2055 + 25 K and lo- 
Strong water-vapor bands have been detected in p Cep cated 0.33 R* above the photosphere. The model proved to con- 
and Betelgeuse by |Danielson, Woolf & Gaustad (1965)| but sistently explain the FLUOR/IOTA K-band measurements and, 
no explanation of their origin was available at that time, after removing the contribution of dust, those of the ISI interfer- 



Jennings & Sada (1998)|detected pure rotational lines of H 2 at ometer at 11.15 pm of [Werner et al. (200"0)| The same conclu- 



12.3 pm in Betelgeuese and Antares and suggested that the water sion was reached on p Cep based on narrow band measurements 
lines may be formed in their atmosphere. Tsuji (2000a, 2000b) in the K band to isolate molecular bands from the continuum. 
The shell has been found to have a 1.32 R+ radius (Perrin et al. 



Send offprint requests to: G. Perrin, e-mail: guy .perrin@obspm. fr 2005). This analysis was also successfully applied to Mira stars 

* Based on observations collected at ESO with the VLTI/MIDI in- (Mennesson et al. 2002, Perrin et al. 2004b) and led to a con- 

strument. sistent understanding of seemingly contradictory measurements 
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Other recent modeling studies of supergiants (Ohnaka 
2004a) and Mira stars (Ohnaka 2004b, Weiner 2004) have shown 
that the near- and mid-infrared H2O spectra and the apparent 
near- and mid-infrared angular sizes of both stellar types can be 
consistently understood in terms of a star with an envelope of 
warm H2O vapor. These results altogether confirm and vindicate 
Tsuji's MOLsphere. In a more recent work, Tsuji (2006) shows 
the consistency of interferometry and spectroscopy-based con- 
clusions. |Verhoelst et al. (2006)] have made a combined analy- 
sis of the near and mid-infrared undispersed interferometric and 
spectroscopic data and showed in addition to H2O, OH and CO 
the likelihood of amorphous alumina in the MOLsphere, a dust 
species with a high sublimation temperature. Ryde et al. (2006) 
proposes an explanation for mid-infrared TEXES spectroscopic 
data which does not require a MOLsphere. Absorption fea- 
tures of water vapor and OH can be reproduced with a cooler 
outer photospheric structure of temperature 3250 K. However, 
this is not compatible with near-infrared spectra according to 



Tsuji (2006) 



Table 1. Reference sources. Diameters are based on photo- 



metric estimates by Dycketal. (1996) with errors revised by 
|Perrin et al. (1998)| Although not directly used for calibration, 
/3 Lep and rj Eri have been used for data selection. 



HD number 


Source 


Spectral type 


Uniform disk 




name 




diameter (mas) 


HD 39400 


56 Ori 


K1.5 lib 


2.25 + 0.12 


HD 49161 


17 Mon 


K4III 


2.65 ±0.13 


HD 36079 


/3Lep 


G5 II 


3.04 ±0.15 


HD 18322 


77 Eri 


Kl III 


2.57 ±0.13 



Although excellent candidates are known as potential con- 
tributors to the opacity of the MOLsphere, other sources may 
be considered such as electron-hydrogen collisions as proposed 
by |Tatebe & Townes (2006 ) for the shell around Mira stars. The 
combination of spectroscopy and interferometry will help in this 
research. In the present paper, spectrally dispersed fringes in the 
N band obtained with MIDI at VLTI allow to evaluate the com- 
pact spherical shell model for a new wavelength range. We first 
present the observations in Sect.|2]and discuss the data reduction. 
Visibility data are fitted with a simple uniform disk + extended 
resolved source model in Sect.[3]to yield normalized visibilities 
free of any extended component. The compact source visibilities 
are then compared to a star + molecular spherical shell model 
in Sect. [4] The shell optical depths are compared to a spectro- 
scopic model in Sect. [5] to derive column densities and identify 
molecules from their signatures. Lastly, the results and their con- 
sequences are discussed in Sect. [6] 

2. Observations and data reduction 

Observations took place at Cerro Paranal in Chile and made use 
of the mid-infrared beamcombiner MIDI (Leinert et al. 2003) of 
the VLTI on the nights of November 8 and 10, 2003. These ob- 
servations were part of a Science Demonstration Time run and 
made use of the UT2-UT3 baseline whose 46.63 m length al- 
lows to sample the first zero region of the visibility function of 
Betelgeuse. 

The MIDI data comprise different sets, two of which are im- 
portant for data reduction. The chopping set provides the aver- 
age flux of the source from which the strong background con- 



tribution has been removed. The fringe set is a collection of 
about a hundred 50 fim long optical path difference scans. At 
the time these observations took place, the MIDI beams were 
not spatially filtered. A 0.6 arcsec slit was used corresponding to 
a 2.5A/D scale at the central wavelength of 10 fim. The beams - 
the two complementary interferometric outputs of the beamcom- 
biner - were dispersed by a prism at a spectral resolution of 30 in 
the direction perpendicular to the slit. No active system was used 
to correct wavefronts except for the tip-tilt correction system to 
maintain the pointing quality. Given that pupils are larger than 
the coherence patch of atmospheric turbulent phase (about 6 m 
for median seeing at Paranal at 10 /urn) and that the beams were 
not filtered, coherence losses due to turbulence fluctuate from 
one scan to another and may prove difficult to calibrate. This is 
the major source of error in measuring visibilities as is discussed 
below. 

An interferogramis the sum of three signals: the background, 
the flux of the source and the fringe modulation. Because the 
beams are not filtered and scintillation can be neglected at 10 /jm 
with an 8.2 m telescope, the flux of the source is constant 
throughout a batch and can be accurately measured with the 
chopping sequence. The quick and large fluctuations of the back- 
ground prevent from simultaneously measuring the source flux 
on interferograms with the MIDI set-up. As a consequence, the 
fringe modulation has to be normalized with the estimate of 
the flux of the source derived from the chopping sequences. To 
do so, the two complementary interferometric outputs are sub- 
tracted to eliminate the background and the instantaneous source 
flux leaving a zero mean signal. The power spectral density 
of normalized interferograms is integrated to yield the squared 
coherence factor of each interferogram in a batch. The power 
spectral density of the sum of noises (background photon noise, 
source photon noise and detector read-out noise) has been sub- 
tracted before the integration to remove the bias on visibility 
modulus due to additive noises. The average of all squared co- 
herence factors in a batch yields an estimate of the squared co- 
herence factor for the batch and their dispersion yields an es- 
timate of the attached uncertainty. A first selection of data oc- 
curs here. Data are rejected when squared coherence factor his- 
tograms strongly depart from a gaussian distribution as this is 
the signature of quick evolution of turbulence statistics during a 
batch thus strongly biasing the average value estimate. 

This procedure is applied to both Betelgeuse and calibrator 
star data. Calibrator star squared coherence factors are divided 
by expected squared visibilities estimated from their uniform 
disk diameter value to compute the squared transfer function. 
When observations of Betelgeuse are bracketted by two calibra- 
tor observations, transfer functions are interpolated at the time 
Betelgeuse was observed as explained in |Perrin (2003)| and er- 
ror bars are computed accordingly. The final squared visibility 
estimate of the source is the ratio of the squared coherence fac- 
tor and of the transfer function and error bars are derived from 
the errors on these two quantities. This procedure is applied to 
each individual spectral channel. It is simpler than what is usu- 
ally done in the near-infrared where turbulence fluctuations are 
filtered with single-mode fibers. Non stationary errors due to tur- 
bulence are impossible to accurately calibrate here and there is 
no need for a sophisticated propagation of noises and uncertain- 
ties. 

The calibration procedure is applied after data have been 
checked for quality. Transfer functions are plotted as a function 
of wavelength for each calibrator. Ideally they should all overlap 
to within noise. Transfer functions are rejected when they clearly 
depart from the average trend. This can be easily detected in as 
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Table 2. Log of November 8, 2003 observations. 
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UT time 


Projected baseline 


Azimuth 1 ' 


opd offset 


Calibrator 1 


UT time 


Calibrator 2 


UT time 




(m) 


o 


fjum) 










06:11:07.00 


37.320050 


38.794453 


250 


56 0ri 


05:29:14.43 


56 Ori 


06:40:55.00 


06:12:07.28 


37.396414 


38.907380 


185 


56 0ri 


05:29:14.43 


56 Ori 


06:40:55.00 


06:16:01.20 


37.691753 


39.335562 


185 


56 Ori 


05:29:14.43 


56 Ori 


06:40:55.00 


07:10:06.99 


41.509800 


43.801921 


180 


56 0ri 


06:42:11.45 


17 Mon 


07:43:28.00 


07:10:59.20 


41.565253 


43.853488 


180 


56 Ori 


06:42:11.45 


17Mon 


07:43:28.00 


08:09:18.00 


44.673489 


46.091170 


170 


17Mon 


07:44:44.45 






08:10:10.20 


44.709776 


46.107820 


170 


17Mon 


07:44:44.45 







counted positive from East towards North 



much as a sudden increase of turbulence or a bad image over- 
lap will cause the transfer function to drop and to drop all the 
more as the wavelength is shorter. Once the selection of trans- 
fer functions is performed, Betelgeuse data are calibrated and 
the same procedure applies to visibility curves as a function of 
wavelength. The application of this procedure has led to reject 
all November 10, 2003 data and part of November 8, 2003. 

All the scans of the first night on Betelgeuse were recorded 
off the zero optical path difference. Fortunately, in dispersed 
mode, the spectral resolution is high enough that the offset is 
within the coherence length of the radiation. Most offsets were 
in the range 170-180 /mi with one batch with an offset of 250 /urn 
(see Table |2j. In the worst case, the offset caused a contrast re- 
duction of 86% at 8 fim and 14% at 13 fim. But for the other 
batches, the contrast reduction was 60% at 8 fim and 7% at 
13 fim. In order to calibrate these losses, we simulated the in- 
terferograms recorded by each pixel of the MIDI detector taking 
into account the dispersion law of the prism and the effects of 
diffraction in the MIDI spectrograph. Losses only depend on the 
characteristics of the instrument and can therefore be reliably 
computed. Laws of contrast loss as a function of wavelength 
were derived for each optical path difference offset and used to 
calibrate all acquisitions. All the calibrator scans were recorded 
at zero optical path difference and no correction was necessary 
for calibrators. 

The entire data reduction has made use of the software de- 
veloped by Paris Observatory for MID0 The list of calibrators 
and the log of calibrated data are given in Table [TJ and Table [2] 
Three different sets of projected baselines have been observed. 
We have averaged data to get three bins corresponding to the 
sets. Averaging has been performed using the statistical errors 
provided by the data reduction software. Each average is a x 2 
estimator: 



x 2 = 



1 



N ■ 



N-l 



Vf) 



2\2 



(T- 



(1) 



The statistical variance of the average visibility is: 



cr 2 = 



N-l 



!=0 



(2) 



The x 2 can be used as an indicator of the correctness of error 
bars and is useful to improve them. If the x 2 is smaller than 1 



1 Software is available through the JMMC website: 
http://mariotti.ujf-grenoble.fr/ 



then statistical error bars are declared correct and or. = a 2 . 

bin 

Otherwise, the error bar on the averaged visibility is corrected 
applying: 



2 2 2 



(3) 



which is equivalent to forcing the x 2 at minimum to 1 This pro- 
cedure is meant to account for non stationary fluctuations of tur- 
bulent errors which are not detected by the data reduction soft- 
ware. This does not perfectly remove biases as these are not 
modeled but at least error bars are increased to better model 
fluctuations. Typical error bars are a fraction of a percent. The 
validity of this approach is assessed in the next section. Given 
the normalization we had to apply to the data, it is difficult to es- 
timate the validity of our error bars as absolute errors. However, 
the derived error bars are internally consistent. The log of aver- 
aged observations is given in Table [3] 

3. Diameter fit 

Betelgeuse is known to be surrounded by a dust envelope whose 
characteristic size of l"(Danchi et al. 1994) is far larger than the 
interfering field of MIDI and decreases the spatial coherence of 
the beams. To first approximation, the source can be modeled 
with an extended and a compact component. The extended com- 
ponent leads to a decrease of the compact component visibility 
by a factor equal to its fractional flux. The study reported here 
only addresses the compact object. We therefore need to normal- 
ize visibilities to clean the compact source from the disturbing 
extended component. 

Visibilities in all spectral channels have been fitted by the 
scaled uniform disk function: 



V 2 (0 UD (A),A,A;Bi)=A : 



2Ji (n0 UD ^) 



7101 



(4) 



where A is the fractional flux of the compact component, 0ud{X) 
its diameter and B, the projected baseline. We assume here that 
the visibility of the extended component is equal to 0. Visibilities 
have been taken close to the first null of the compact source visi- 
bility function. The effect of finite bandwidth has to be taken into 
account as the estimated squared visibilities do not go through 
0. For comparison with observations, the above model has there- 
fore been integrated in each spectral channel as described in 
|Perrin et al. (2004a)| 

75% of the fits have^ 2 values smaller than 1 showing that the 
calibration and the correction of error bars are correct in average. 
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Table 3. November 8, 2003 averaged observations. 



Mean UT time 


Mean projected baseline 


Mean azimuth" 


Calibrator 1 


UT time 


Calibrator 2 


UT time 




(m) 


O 










06:13:36.12 


37.508548 


39.906969 


56 0ri 


05:29:14.43 


56 Ori 


06:40:55.00 


07:10:33.09 


41.537526 


43.827705 


56 Ori 


06:42:11.45 


17Mon 


07:43:28.00 


08:09:44.10 


44.691632 


46.099496 


17Mon 


07:44:44.45 







counted positive from East towards North 
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Fig. 1. Examples of MIDI visibility fits by the scaled uniform disk visibility function (solid line). Fits with the spectroscopic model 
of Sect. [53] with a 1520 K layer are also presented (dashed line). 



Examples of fits are presented in Fig. [TJ Data below 9.18 fim 
sample the second lobe of the visibility function of the compact 
component whereas data above 9.56 fim sample the first lobe. 
Confidence in the second lobe arises from the slope of the visi- 
bilities which at all wavelengths is positive. Because the source 
is probably limb-darkened, the fractional flux is likely to be un- 
derestimated in the second lobe as a reduced slope due to limb- 
darkening will be compensated by a lower A with the uniform 
disk model. Visibilities between 9.18 and 9.56 /urn are difficult 
to fit as they correspond to the first zero of the visibility func- 
tion and are very sensitive to errors. They are not used in the 
following. 



Fig. [2] shows the fractional flux of the compact component 
and its apparent diameter as a function of wavelength. The graph 
can be separated into three regions. Below 9.18 /urn, the frac- 
tional flux of the compact component is very large. However, 
because of the bias due to the unknown limb darkening, the frac- 
tional flux in this range is probably underestimated and makes its 
interpretation difficult. A trough is clearly visible around 9.7 //m 
which could be the signature of silicates in the bright resolved 
dusty environment. A smaller trough is also present between 1 1 
and 12 pm which could be due to alumina (AFO3) - as it may 
be present in low-density shells (Sedlmayr & Kriiger 1997) - or 
other species. The fractional flux above 10 pm is roughly con- 
stant but significantly smaller than 1 which could possibly be the 
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Fig. 2. Fit of the MIDI data by the scaled uniform disk visibility function. Left is the fractional flux of the star A. Right is the uniform 
disk diameter 0ud- The MIDI diameters are compared to the ISI diameters at 1 1.086, 1 1.149 and 11.171 fim (Weiner et al. 2003). 
The 1 1 . 149 /im ISI value is the average of measurements taken at different epochs and the error bar is the standard deviation of these 
measurements. 



onset of the signature of the 18/jm silicate feature. The apparent 
size below 9.56 fim is comparable to the size measured in the 
visible probably indicating that the same regions are sampled in 
these respective wavelength ranges. The optical depth would be 
large in the mid-infrared and the radiation could be scattered in 
the visible. Above 9.56 fim the apparent size steadily increases. 

The complex behavior of the apparent size of the star can 
be most simply understood as a result of wavelength dependent 
opacity in an extended atmosphere. 

4. Spherical shell model 

We have shown in the previous section that the apparent diame- 
ter of the compact source, the star embedded in the dusty shell, 
varies with wavelength. In |Perrin et al. (2004a)| we had proposed 
a model to consistently explain the star diameter difference be- 
tween the near-infrared (K and L bands) and the mid-infrared 
(11.15 fim). The model consisted of a photosphere surrounded 
by a spherical gaseous shell possibly containing H2O and SiO. 
This fairly simple model is consistent with the more elaborated 
model of the MOLsphere of Tsuji (2006)| The independent set 
of data presented in this paper offers an excellent opportunity to 
test and refine this model. 

The MIDI visibilities have been rescaled by the reciprocal of 
the fractional flux factor A of Eq. [4] to cancel the visibility loss 
due to the more extended environment. All wavelengths have 
then been fitted by the star + spherical shell model whose spatial 
surface brightness distribution is defined by: 



(5) 



for// > yi -(0*/0 sh eii) 2 and: 
/(^//) = B(^,r shell )[l-e- 2T ^] 



(6) 



otherwise. 

fi = cos 6 with 6 the angle between the radius vector and the 
direction to the observer. The parameters of the model are the 



diameter of the photosphere 0* and its temperature T*, the di- 
ameter and temperature of the shell s h e ii and r s h e ii, and the shell 
optical depth t a , the latter being the only wavelength-dependent 
parameter of the model. 

We have considered that the photosphere diameter value has 
been well established from near-infrared observations. This has 
allowed us to fix two parameters of the model: the diameter of 
the photosphere and the temperature. We have taken values from 
|Perrin et al. (2004a)| 0* =43.71 mas and = 3641 K. The 
diameter is the average of the two limb-darkened disk diameter 
estimates of that paper. A consistent limb-darkened diameter 
has been measured in the H band with visibilities in the first and 
second lobes ( |Haubois et al. (2006)[ ). In the following, we use 
this near-infrared limb-darkened diameter as an N-band uniform 
disk diameter as limb-darkening in this latter wavelength range 
is supposedly much smaller. This is however a limitation as we 
do not measure any limb-darkening in the N band and interpret 
apparent limb-darkening as due to the opacity of the shell. As 
is discussed later, this is mostly a concern in the lower part 
of the wavelength range where data are only available in the 
second lobe of the visibility function. Fixing the photosphere 
temperature and diameter parameters was anyway necessary as 
all parameters may be well constrained by the curvature of the 
visibility function only if the first lobe of the visibility function 
is well sampled with both high and low visibilities which is not 
the case with this set of data. Our fixing of these two major 
parameters at these particular values is however justified as 
|Perrin et al. (2004a)| have shown that they are in agreement with 
the bolometric flux and spectral type of the object. The si- 
multaneous fit of all wavelengths provides the following results: 

0sheii = 62.30 + 0.13 mas 
r shell = 1598 + 20 K 
X 1 = 6.29 

In this fit, since the temperature and the size of the photosphere 
are fixed, the shell temperature reaches a single minimum at op- 
timum value and does not require the total flux of the observed 
object to be determined in order to mitigate the issue of degener- 
acy as is the case when the parameters of the photosphere are to 
be determined. However, a better solution is found if two sepa- 
rate wavelength domains are fitted independantly corresponding 




Fig. 3. Result of the fit of the normalized visibilities by the photosphere + shell model. The results are obtained with two independent 
fits shortward of 9. 1 8 /mi and longward of 9.56 /mi. The left panel is the optical depth of the shells as a function of wavelength. The 
model spectrum is compared to the measured and normalized spectrum (to account for the dust around the star, see text) in the right 
panel. The shell temperatures are adjusted to better fit the measured spectrum. 



to two different altitudes for the shell: shortward of 9. 18 /mi and 
longward of 9.56 //m as could be guessed from Fig. [2] But in 
this case the shell temperatures are not well constrained and it is 
necessary to force the model spectrum to be consistent with the 
measured spectrum to get accurate and realistic temperatures. 
The following optimum parameters have been found for the two 
domains: 

A< 9.18 /mi A > 9.56pm 

she ii = 62.50 ± 0.50 mas sheU = 57.25 + 0.03 mas 

r shell = 1570 + 150 K r shell = 1520 + 400 K 

with a global reduced^ 2 of x 2 = 3.24. 

The temperatures found for the two shells are very consistent 
with the temperature found with the single-shell model. Because 
the uncertainties are large, it is not possible to discuss the relative 
values of the shell temperatures. We only consider the two-shell 
model in the following. 

The optical depth is plotted as a function of wavelength in the 
left panel of Fig. [3] Some values are poorly defined in both the 
lower and upper wavelength sub-ranges. Below 9.18 /mi, data 
have been obtained in the second visibility lobe and this spa- 
tial frequency range is very sensitive to limb darkening. In this 
case the fractional flux may have been underestimated and has 
led to large visibilities, sometimes higher than the summit of the 
second lobe. When this happens, the optical depth and the error 
bar have to be large. Limb darkening in the second visibility 
lobe also has an impact on the derived diameter. The compari- 
son of the diameter found for the shell below 9.18 /mi and the 
maximum apparent diameter found in Section [3] may be found 
surprising as the shell diameter is smaller by 2 mas. This is most 
probably due to a bias in the fit of the second visibility lobe data 
by a uniform disk function. The size of the object is very sensi- 
tive - to this level of accuracy - to the assumed amount of limb 
darkening. With the photosphere + shell model, the photosphere 
is naturally limb darkened when the shell has a non zero optical 
depth. Visibilities in the first lobe would certainly allow to de- 
rive a better estimate of the apparent diameter more in agreement 
with the derived shell diameter. Data taken in the first visibility 
lobe above 9.56 /mi are far less sensitive to this effect. However 



a few values are large with large error bars. This can be at- 
tributed to calibration issues in this upper part of the wavelength 
range where sensitivity to background is much larger. These val- 
ues are not to be taken into account for the interpretation of the 
optical depth. 

The apparent diameter found at 11.15 /mi is quite consis- 
tent with that measured with ISI by |Weiner et al. (2000)] re- 
spectively 53.75 + 0.32 mas versus 54.74 + 0.34 mas. The dif- 
ference may arise from different reasons one of which is that 
the bandwidth of ISI is by far much narrower than MIDI'S, re- 
spectively about 2 nm versus 100 nm. The two shells we find 
are just above the shell of |Perrin et al. ( 2004a) whose diameter 
was 55.78 mas. These altitudes are satisfyingly consistent as the 
data sets are difficult to compare and could be improved in ei- 
ther case. There is also a good agreement between these shells 
and the (62.35 mas, 1750 K) shell of |Verhoelst et al. (2006)| 

The spectrum of Betelgeuse measured by MIDI is also pre- 
sented in Fig. [3] More exactely, we have multiplied the spec- 
trum measured by [Verhoelst et al. (2006)| with the same set of 
MIDI data by the fractional flux of the central object derived 
in Section [3] and presented in Fig. [2] in order to eliminate the 
contribution of dust to the spectrum. The error bars have been 
estimated to 0.01 x 10~ 13 Wcm^m -1 to account for the uncer- 
tainty on the flux ratio. Error bars for the temperatures of the 
two shells have been determined by matching the model flux un- 
certainty to that of the measured spectrum. On the same figure, 
upper and lower limits for the optical depth are plotted (dashed 
lines). They correpond to the optical depth derived for ±l<x val- 
ues of the temperature as constrained by the spectrum. They give 
an idea of the probable range of values to be expected for the op- 
tical depth as this range is larger than the purely statistical error 
bar derived from the x 1 analysis. 



5. Spectroscopic modeling 

In this section, we derive the composition of the layer(s) found in 
the previous section. As discussed by Scholz (2001), low spec- 
tral resolution visibilities of extended atmospheres with strong 
molecular opacity cannot be correctly modelled in any other 
way than by computation of the model and the corresponding 
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visibilities at high resolution (resolving the individual lines), be- 
fore convolving those resulting visibilities with the instrumental 
spectral response. Yet, our data set comprises too few baseline 
coverage and sometimes in a very sensitive spatial frequency 
range (the second lobe of the visibility function below 9.56 fim) 
to perform this direct and rigorous analysis. We therefore first 
do a comparison of molecular and dust opacities to the empiri- 
cally derived opacity profile of Sect. |4] over the N-band, the ge- 
ometrical model being used as a proxy of the star. This allows 
the identification of possible constituents and to estimate col- 
umn densities. We also perform a direct confrontation between 
model visibilities and observed visibilities in Sect. l5.3l and study 
the limits of this approach in the context of our current dataset. 
In the following, we assume a single layer with a temperature of 
1520 K and reaching from 1.32 to 1.42 R*. This is a more realis- 
tic representation of the two infinitesimally thin layers found in 
the previous section and allows the use of a homogeneous dis- 
tribution of the different components^ We have chosen to derive 
error bars on the column density based on the observational er- 
ror in the optical depth determination under the assumption of a 
single-temperature layer. Assuming a single temperature across 
the N band yielded an uncertainty on this temperature of only 
20 K in the previous section whose impact on t uncertainties is 
negligible. We did however rescale the error bars on the t mea- 
surements uniformly in order to obtain a x 2 r — 1 for the best-fit 
model. 

5.1. Possible sources of opacity 

From the molecules commonly found in CSE's, H2O and SiO 
are known to have strong opacity bands in the MIDI wavelength 
range (see e.g. Decin 2000). SiO is an absorber below 10 fim 
whereas H2O absorbs across the entire range. Other molecules 
such as silane and ammonia may provide high opacities in 
carbon-rich stars and supergiants (see Monnier et al. 2000 and 
references therein) but were not found to be good candidates for 
our analysis. In this section, we compare the observed opacity 
profile (Fig. with optical depth profiles which we computed 
for H2O and SiO at the temperature derived in Sect. [4] 
For H2O, we used the nasa ames linelist and partition func- 
tion (Partridge & Schwenkel997), restricted to lines with 
l°g ^/ ~ X® ^ ~9, with g the ground level statistical weight, 
/ the oscillator strength, 9 = 5040/3500 and x tne excitation 
energy in eV. For a justification of the use of this particular 
linelist and of the applied cut-off in line strength, we refer to 
Cami (2002), Van Malderen (2003) and Decin (2000). For SiO, 
we used the line list of Langhoff & Bauschlicher (1993), which 
is the most complete list to date (Decin 2000). The polynomial 
expansion of the appropriate partition function was taken from 
Sauval & Tatum (1984). Although very high rotational and vi- 
bration levels are included in the line list, SiO causes opacity 
in the N band mostly in its fundamental mode at 8.12 /im, and 
in two hot resonance bands at 8.20 and 8.28 /mi. The micro- 
turbulent velocity was taken to be 3 kms _I , a typical value for 
molecules in AGB stars (Aringer et al. 2002) which is probably 
also valid for late-type supergiants. 

The only dust species which can survive at these high temper- 
atures and which provides significant opacity at these wave- 
lengths is alumina (AI2O3). We used the optical constants mea- 
sured by Begemann et al. (1997) and Koike et al. (1995). A grain 



2 A comparison with a model containing 2 separate layers does not 
show any significant differences at the spectral resolution of MIDI and 
ISO-SWS. 
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Fig. 4. Optical depths observed with MIDI (black diamonds) 
and high-resolution opacity profiles of water (black), SiO (grey) 
and alumina (solid grey line). 

Table 4. Derived column densities of H2O, SiO and alumina. 



Component 


N [cm" 2 ] 


H 2 


7.1 ±4.7 x 10 19 


SiO 


4.0 ± 1.1 x 10 20 


AI2O3 a 


2.4 + 0.5 x 10 15 



a grain size of 0.005 yum and density of 4gcirT 3 were assumed 
with opacities computed in the Mie approximation (spherical 
grains) 



size of 0.005 fim and a density of 4 gem -3 were assumed with 
opacities computed in the Mie approximation (spherical grains). 
In Fig. [4] we show the observed opacities across the N band, and 
the high-resolution opacity profiles of H2O, SiO and alumina. 
We see that the high opacities below 9.18/im suggest the pres- 
ence of a significant column of SiO. Beyond 9.56 /mi, the in- 
crease of opacity with wavelength is compatible with the pres- 
ence of both H2O and alumina. 

5.2. Modelling the observed opacity profile 

We attempt a first estimation of the layer composition by direct 
comparison of opacities. 

The computed opacity profiles were rebinned to the resolution 
of the MIDI observations (R~30) and column densities were de- 
termined by searching for the best agreement with the observed 
optical depths. The derived column densities are presented in 
Table H] and the result is shown in Fig. [5] To determine realistic 
uncertainties on the derived column densities, we rescaled the 
errors on the observed optical depths so as to obtain ax 2 ~ F 

We find that the observed optical depths can be reasonably 
reproduced with H2O and SiO column densities of the order of 
10 20 cirT 2 and an alumina column density of 5 mgcrrT 2 . The 
near- to mid-IR implications of such a layer, which can be tested 
against the ISO-SWS spectrum, are shown in Fig. [6] In general, 
our model represents the ISO-SWS spectrum reasonably well. 
The overestimation of the absorption by H2O at 2.5 //m suggests 
that our layer temperature may be about 50 K too low. Still, this 
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Fig. 5. Observed optical depths (black diamonds: optimum tem- 
peratures, grey diamonds: upper and lower limits on tempera- 
tures) compared to the best fitting model opacity profile (thick 
grey solid curve). 



Fig. 7. Observed and modeled MIDI visibilities are shown in 
red, green and blue for the 37.5 m, 41.5 m and 44.7 m projected 
baselines respectively. The black crosses show the best-fit scal- 
ing factor for each wavelength. 
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Fig. 6. Upper panel: the near- to mid-IR ISO-SWS spectrum 
(black) and the model with a molecular and dusty (H2O, SiO and 
alumina) layer (grey). The dusty, detached shell fully resolved 
by MIDI is not modeled here hence the discrepancy at the sili- 
cate emission peak near 9.7 jum. Bottom panel: opacity profiles 
of H2O, SiO and alumina (dotted line) for the temperatures and 
column densities in our model. 



model is in agreement with that proposed by Verhoelst et al. 
(2006). 

5.3. Confrontation with the visibilities 

Since our approach in the previous section is not rigorously cor- 
rect, we test our model against the observed visibilities. 

5.3.1 . With the column densities of Sect. |5.2| and the 
geometrical model parameters 

We first computed visibilities across the N-band at the MIDI res- 
olution based on the column densities derived from the t fitting 
in the previous section and compare these to the observations in 
Fig. [7] We allowed for a new determination of the scaling factor 



due to the flux of the resolved component that takes into account 
the limb-darkening of the photosphere by the molecular shell. 
We see that the match in the first lobe of the visibility curve, 
i.e. the wavelengths longward of 9.56 //m, is quite good. The 
agreement is much poorer at shorter wavelengths. This is due 
to insufficient opacity of the layer at these wavelengths: the SiO 
lines (at this layer temperature) are too sparsely spaced so the 
visibility is still dominated by photospheric photons, an effect 
which can't be overcome by increasing the column density. This 
is much less a problem for H2O which has many more closely 
spaced lines. Moreover the opacity of our model at longer wave- 
lenghts is dominated by the continuum opacity of alumina. 
The scaling function clearly shows the silicate features at 9.7 
and 18 fim in absorption, indicating that the resolved component 
is, as expected, the olivine-rich circumstellar dust shell. This is 
also compatible with absorption by alumina or other species at 
11.3 fim. 

5.3.2. With column densties and layer diameter as free 
parameters 

To improve on our determination of the layer composition, we 
searched a grid of models for the best agreement with the ob- 
served visibilities. Our grid ranges from 57 to 62 mas in layer 
diameter (1.30 to 1.42 R* in layer radius) with 1 mas step- 
size, 1 x 10 19 to 1 x 10 22 cirT 2 in H2O column density with a 
1 x 10 19 ' 20 ' 21 stepsize, 1 x 10 20 to 1 x 10 24 cnT 2 in SiO column 
density and 0.05 to 8 mg cirT 2 in alumina column density. The 
lowest^ 2 is found for a layer of 58 mas in diameter (1.33 R* in 
radius), with column densities of 3 x 10 21 crcT 2 , 1 x 10 20 cm~ 2 
and 0. 1 mg crrT 2 respectively. This model is shown in Fig. [8] The 
X 2 r (rescaled to 1 for the best model) as a function of the different 
parameters is shown in Fig. [9] We can see that the water col- 
umn density is reasonably well constrained. The SiO has in fact 
no noticable effect on the low-resolution visibilities. Alumina 
can be present, but we find an upper limit of approximately 1 mg 
cm 2 . In contrast to the empirical modelling presented in Sect. [4] 
we find here a very poor constraint on the layer diameter. 
Since a water column density of the order of 10 21 cnT 2 is incom- 
patible with both near-IR - interferometric and spectroscopic - 
observations (too much opacity) and with the high-resolution 
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Fig. 8. Best matching model in our grid of 1520 K layers. Fig. 10. Model with a layer at 2000 K and containing H2O, SiO 
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Fig. 9. x 2 r (rescaled to 1 for the best model) as a function of the 
different parameters. 



ISI observation at 1 1 /mi (too little opacity), we may consider 
a higher layer temperature, e.g. 2000K: this strengthens the hot 
resonances of SiO, resulting in lines which are spaced much 
closer. The increased SiO opacity below 9.18 /mi could possi- 
bly reduce the required amount of H2O and increase the required 
column density of alumina. 

A fit to the observations of a model with a layer at 2000 K 
containing H2O, SiO and ALO3 is shown in Fig. \TU\ We find 
the best agreement with column densities of 7.4 x 10 20 cm~ 2 , 
8. Ox 10 22 crrT 2 and 4.2 x 10~ 6 g crrT 2 respectively and a layer di- 
ameter of 59.13 mas (1.35/?* in radius). Indeed, we find that this 
does increase the SiO column density but lowers H2O column 
density not enough and the alumina column too much. It there- 
fore brings no improvement in the near-IR and high-resolution 
1 1 /mi properties of our model. 

At this stage we conclude that a rigorous modelling of our data 
is not possible and that at least a better sampling of the visibility 
function at smaller spatial frequencies shortward of 9. 18 /mi is 
required. 



6. Discussion 

6.1. The MOLsphere 

These are the very first interferometric data that allow to study 
the characteristics of Betelgeuse over a wide mid-infrared 
spectral range. The only other mid-infrared data available so 
far have been collected by the ISI interferometer. Comparison 
of our derived diameters with that of ISI at 11.086, 11.149 and 
11.171 /mi (Weiner et al. 2003) is excellent as Figure [2] shows 
although the ISI spectral resolution is an order of magnitude 
higher than ours. The strength of the MIDI data is their spectral 
coverage which allow us to go further in the spatial and spectral 
understanding of the source. We derive the optical depth 
of the environment above the photosphere as a function of 
wavelength. Absorption is clearly not due to a continuum source 
of opacity across the N band as the optical depth drops above 
9 fim. Several types of absorbers are required to reproduce the 
measured optical depth. 

The ad-hoc model of the MOLsphere is well adapted to capture 
most of the source characteristics both interferometric ally 
and spectroscopically. Potential constituents are identified. 
Spectro-interferometric data are well modeled with H2O, AI2O3 
and SiO in the MOLsphere. This confirms the presence of 
dust species which can survive close to the star thanks to 
their glassy nature - the low opacity of alumina at optical 
and ultraviolet wavelengths prevents strong heating by the 
star radiation (Woitke 2006). Furthermore, the presence of 
gaseous phase SiO in the MOLsphere strengthens the scenario 
of |Verhoelst et al. (2006)| for dust formation, whereby alumina 
may condense first and act as a nucleation site for silicate dust 
which can condense further from the star where the temperature 
drops below 1000 K. The current interferometric data set is 
not enough to allow a more complex model and we had to 
assume a geometrically thin shell around the star. In particular 
as underlined in Section[3]a more complete set of data is needed 
both at shorter and somewhat longer baselines in order to better 
measure limb darkening and better measure the optical depth 
of the shells above the photosphere. This would allow to test 
this model further and possibly discuss more complex spherical 
structures and temperature-density radial variations. 
With this rich data set, it is interesting to review our results 
in the light of other available studies of the MOLsphere. Our 
modeling mostly accounts for the ISO-SWS spectrum at a low 
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resolution level. However, even at low resolution the agreement 
is not perfect (for example between 2.4 and 3.2 /mi) and the 
too simple description of the MOLsphere in our model is 
potentially an issue. Apart from this discrepancy, the agreement 
is globally very good from the upper near-infrared range up 
to the mid-infrared (the silicate dust infrared excess is not 
accounted for in our modeling since the dusty environment has 
been completely resolved out). 

|Ohnaka (2004a)| has modeled a MOLsphere containing water 
vapor only. Similar parameters to ours were used to model 
the star (fl* = 650 R and T eff = 3600 K) but the MOLsphere 
description is a little different: it extends from the photosphere 
up to R mo i = IA5R+. A temperature r mo i = 2050 K and a 
water vapor column density of 3 ± 1 x 10 20 cnT 2 are found. 
In Tsuji (2006) a slightly different description is used. The 
MOLsphere is characterized by its inner radius R[ n above 
which the density is decreasing up to the outer radius R out . 
The same parameters are used to describe the star. Different 
sets of parameters are investigated for the MOLsphere and the 
preferred solution is r mo i = 2250 K with R ln = 1.31/?* and a 
water vapor column density of 1 x 10 20 cirT 2 . The two models 
are very consistent as far as water vapor is concerned. Given 
the different geometries considered by the two authors (R mo i is 
the upper limit of the MOLsphere in one case whereas /?; n is 
the lower limit in the other one), one expects the Ohnaka radius 
to be larger than the Tsuji radius which is exactly as they are. 
The two models define a range of radius for the MOLsphere 
of 1.31 to 1.45/?* which is exactly the range considered in our 
spectroscopic model. Column densities are also very consistent 
to ours. We however find a lower temperature range although, 
as previously stated, the uncertainty on this particular parameter 
is somewhat difficult to define with the current data set. Despite 
this difference, we conclude the different modelings are quite 
consistent and we bring evidence for new constituents of the 
MOLsphere. This is a strong support for the MOLsphere model 



and we agree with Tsuji (2006) that the photospheric origin 
of all water vapor lines suggested by Ryde et al. (2006) is not 



consistent with both spectroscopic and interferometric data. 

It is interesting to compare the results of applying the star- 
shell model to different stellar types. The earlier near-infrared 
study of Mira stars (Perrin et al. 2004b) determined model pa- 
rameters which differ significantly in several respects from the 
model of a Ori. The Mira photospheres are cooler by ~400 K. 
The Mira shell radii are twice the stellar radius (compared to 
~1.4 for a Ori), and yet the Mira shell model temperatures 
are similar to or slightly higher than the Betelgeuse shell 
temperature. These differences highlight the fact that, although 
a similar model geometry can be applied to the two cases, they 
are physically very distinct systems. 

Obviously, one of the next steps will be to better reproduce 
the spectra at a level for which there is an agreement between 
the model and the spectral features at the resolution of MIDI. 
A better sampling of the visibility function is required for that 
to reduce the number of a priori inputs to the model (such as 
the star photospheric radius and temperature) and to determine 
the intrinsic limb darkening of the photosphere and to allow the 
more rigorous computation of the model. 

6.2. Mass loss 

There are many questions concerning the mass loss process 
in stars, and different processes may be at work in different 



environments. In a star such as a Orionis, a persistent problem 
concerns the location of dust nucleation and initial acceleration 
of mass, since a classical hydrostatic atmosphere fails to offer a 
suitable location for dust formation. 

The new evidence reported here, integrated with other recent 
results, brings together important clues which may guide the 
way to a preliminary understanding of dust formation in the 
atmosphere of a Orionis. The H2O column density we measured 
can be converted to a total gas column density (assuming 
«h,o/«h = 10~ 7 according to Jennings & Sada (1998) , and 
a layer thickness of rougly 10% of its radius) which yields a 
mass density of 10~ 13 gcm~ 3 , which is typical for the wind 
acceleration/dust condensation regions in sophisticated models 
(see e.g. Fig. 1 in Woitke 2006). Such a gas distribution, not 
predicted by static models, is therefore dense enough to support 
dust formation chemistry. The low inferred temperature of the 
molecular shell (~ 1600 K) is in the range of conditions believed 
required to permit grain nucleation and growth for one or 
more species which have particularly favorable thermochemical 
properties. 

Finally, a third piece of the puzzle is the tentative identification 
of emission from a condensed grain material (here identified 
with alumina) which is associated spectroscopically and spa- 
tially with the dense, cool region. 

The star-shell model is proving to be a fruitful tool for a 
preliminary reconnaissance of mass loss systems. Not surpris- 
ingly, the shell parameters found for stars of different types 
(red supergiants, Mira stars) are quite different. As it becomes 
possible to estimate molecular compositions and densities with 
increasing confidence, the results may be a guide to the likely 
differing processes at work. 

Clearly many points in this simple description require study, 
including the energy balance and equlibrium temperature of the 
material, the temperature distribution in the shell, and the evo- 
lution of chemical equilibrium. Many questions remain, such 
as the mechanism for establishing the high density, the dynam- 
ics of dust and gas acceleration, and the role of time variations 
and inhomogeneities. Spectrally resolved mid-IR interferome- 
tery appears well suited to observationally constrain such inves- 
tigations. 
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